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Abstract 

The framework of the Perturbed Static Path Approximation (PSPA) is used to 
calculate the partition function of a finite Fermi system from a Hamiltonian with a 
separable two body interaction. Therein, the collective degree of freedom is introduced 
in self-consistent fashion through a Hubbard-Stratonovich transformation. In this way 
all transport coefficients which dominate the decay of a meta-stable system are defined 
and calculated microscopically. Otherwise the same formalism is applied as in the 
Caldeira-Leggett model to deduce the decay rate from the free energy above the so 
called crossover temperature Tq. 
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1 Introduction 



Commonly quantum versions of the decay rate of damped meta-stable systems are 
treated on the basis of the Caldeira-Leggett model For an exhaustive overview on 
this topic we may refer to the review articles by P. Hanggi et al. Q, G.-L. Ingold in 
or the text books by U. Weiss Q. Unfortunately, this model assumes some simplified 
coupling to a linear heat bath. Furthermore it does not make any predictions about 
the input of such important quantities as the potential energy and the inertia, which 
need to be chosen on an entire phenomenological level. These features do not allow 
the model to be applied to self-bound Fermi systems. There, one would like to see 
the collective variables introduced in some self-consistent fashion, with a microscopic 
treatment of all transport coefficients alike. 

One possible attempt to overcome these deficiencies is formulated in (with ref- 
erences to earlier papers) in connection to nuclear physics. It is based on a quantal 
transport equation which is derived within a locally harmonic approximation exploit- 
ing linear response theory. It is this approximation which allows one to treat a more 
complicated coupling between the collective variable and the intrinsic degrees of free- 
dom. A transport equation necessarily describes evolution in real time. Therefore, in 
barrier regions quantum effects can be accounted for only above a critical temperature 
Tc, which is larger than the so called crossover temperature Tq which one encounters 
for imaginary time propagation B. As one knows, the same feature holds true also 
for the Caldeira-Leggett model Another disadvantage of the derivation of this 

transport equation mentioned is that it bases on the deformed shell model. Surely, 
it allows one to calculate all transport coefficients on the same footing. But as one 
does not start from a genuine two body interaction, self-consistency is handled on a 
semi-microscopic level only. 

It is the aim of the present paper to do first steps to overcome these deficiencies. 
This is possible by adapting a previously developed formalism to evaluate the partition 
function for bound systems with separable two body interactions. One starts from the 
so called Static Path Approximation (SPA) functional integrals as an approximation 
to the classical or high temperature limit |^ |l^, Then small scale fluctuations 
around this static path are treated to second order, such that quantum effects come 
in through local RPA. In the literature this is referred to as RPA-SPA in |12], the 
Perturbed Static Path Approximation (PSPA) in [O] (the name which we are taking 



over) or Correlated Static Path Approximation (CSPA) in [14|. 



2 Partition function of a finite Fermi system 

Finally, we are interested to generalize the formulas of dissipative tunneling to a system 
where the collective degrees of freedom are introduced self-consistently. The simplest 
Hamiltonian which may serve this purpose is of the following structure 

H = Ho + ^FF, (1) 

with (hermitian) one body operators Hq and F. The product FF mimics an effec- 
tive separable two body interaction. For isoscalar modes, the case we have in mind 
predominantly, the coupling constant k is negative [^]. As we shall see later, the F 
describes one collective degree of freedom. The ansatz (||) should be considered to 
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define a microscopic model for just tliis collective mode we want to address to |15]. 
Neglecting spin and isospin degrees of freedom, a general two body interaction may be 
written as a sum of separable terms 

= i?o + ^5^A;ii;F,. (2) 

i 

For instance, one might exploit an expansion into multipole operators. In case the 
latter are not hermitian the product must be replaced by F^Fi (see e.g. section 4.4.7 
of 11). 

2.1 The general form of the partition function 

The partition function of the grand canonical ensemble reads 

Z{I5) = Tr exp {-I5{H - /ii)) =TiU , (3) 



where /3 = l/T is the inverse temperature and H is the Hamiltonian (|^). The chemical 
potential ^jl is needed in order to keep the particle number [A) fixed on average. It 
would be more appropriate to work with truly fixed particle number. But as we are 
mainly interested in the dependence of transport properties on excitation energy or 
temperature in this paper, this simplification should be accepted. An exact evaluation 
of dH) is prohibited by the presence of the two body interaction. Treating the latter 
in mean field approximation facilitates the calculation greatly. A convenient technical 
tool to incorporate this approximation is to use functional integrals (in imaginary time 
propagation) with which fluctuations about the mean field may be treated as well. 
An elegant form of handling this problem is given through the Hubbard-Stratonovich 



transformation |18|, by which the collective variable ^(t) is introduced. To keep the 



present exposition as short as possible we safe ourselves from repeating the derivations 



of - |14| but simply state the basic results which will then serve as the starting point 



for our generalizations. Mind, however, that the notation has been adapted to that 



used in transport theory ||19|. 



After introducing the Fourier expansion of the collective variable 

<i{t) = qo + ^qr exp(if^r) , (4) 

where the so called Matsubara frequencies 
27r 27r ^ 

z^r = 177; r = — rT with r = ±1, it 2, it 3 . . . (5) 

hp h 

(in units with k-Q = 1) have been used, the partition function may be written in the 
following form within the PSPA (see eq.(21) of llJ 



^(/^) = \/ r^dqo e^^o q,) C(/3, q^) . (6) 



oo 



-2TTk 

Here 

z{p, qo) = Tr exp (-/3{ho{qo) - M)) =11(1 + exp (-/3(eKgo) - ^ (7) 



3 



is the grand canonical partition function belonging to the static part of the Hamiltonian 
in mean field approximation 



hoiqo) = Ho + Fqo (8) 

which is simply a sum over one body operators. The corresponding one body Schroding- 
er equation at some given qq reads: 

hoiqo) \l{qo)) = ei{qo)\l{qo)) (9) 

The appearance of the qo reflects the static version of the self consistency relation for 
the mean field, 

q = k{F) , (10) 

which relates the collective variable q to the expectation value of the operator F. The 
exponent in the first factor of the integrand of (^ is easily understood as to represent 
the static part of the correction —k{F)'^/2 which the energy (ho) in the independent 
particle picture gets from the two body interaction. Neglecting the factor C one obtains 
the partition function in Static Path Approximation (SPA) 10, ^] 

On the very right of ( pT] ) the symbol J^^^^(P,qo) has been introduced to represent a 
free energy. It is not the one of the total system (or total nucleus) which would be 
given by the relation 

Z{P) = exp i-mP)) , (12) 

when the partition function is identified as Z{P) = Z^^^{f3). Rather, the !F^^^{P,qo) 
represents the free energy of the system of nucleons whose mean field is kept fixed at 
the go- In 8- common language of transport theory one would call it the free energy of 
the intrinsic degrees of freedom. 

So far any contribution from the dynamics in the collective variable g(r) has been 
neglected. Formally this may be accounted for by writing the correction factor as the 
following path integral |13] 

C(/3,go) = / -D'q exp f f V|grP+ln(^),o ) (13) 



with the measure 

V'q = lim TT — - dReiqr) dlm(g^) (14) 

Ne=hl3 r=l 

In ([l^) there appears the thermal expectation value of an evolution operator lAq which 
can be expressed by the following (imaginary) time-ordered product 



i^<i)qo = Tr V ^ ( exp(-/3/io(go)) '?'exp 



his 

dThi{T,qr) 



(15) 
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The Hamiltonian 



hi{T,qr)=F{T)5q{T) (16) 

may be understood as the time dependent correction to the static mean field given in 
(^). Here, the time-dependence of the operator part is defined as 

F{t) = e'^o^'^'')^/'* F e~^o(''o)'^/^ , (17) 

which means through the interaction picture based on the Hamiltonian hQ{qQ) of (|8|) and 
thus depends on the static qq. The fluctuation of the collective variable 6q{T) = ^(t)— go 
in (16) is related to the fluctuating mean field through (0). It may be noted that the r- 
dependence of the c-number q is meant to be the correct one, not that of any interaction 
picture. 

The partition function (^) may finally be written in the following compact form 

ZiP) = f°°dqo exp(-/3^(/3, qo)) (18) 

if again one uses the concept of the " intrinsic free energy" which now is given by 

nP,Qo) = -^Qo - ^ ln^(/3,go) - ^ lnC(/3,go) ■ (19) 

2.2 The Perturbed Static Path Approximation (PSPA) 

We are now going to evaluate the general formula (|l^) within the so called PSPA. It 
is defined as that approximation in which the exponent appearing in (^) is expanded 
to second order in the qr- This leads to the common Gaussian approximation which is 
known to be related to the semi-classical limit. Following |13] one may write 

l^K)9o' = ^ E 'iris / / d^e--e--(TF(r)F(a)),„ , (20) 

with the r-dependence of the operators as defined in (|l7|). Likewise, according to 
(15), the expectation value is to be calculated with the density operator corresponding 
to the same unperturbed Hamiltonian /io('7o)- It is this feature which will allow us to 
introduce and work with response functions. As we shall see below, this is of advantage 
for at least two reasons, which in a sense are related to each other. The final result, say 
for the decay rate of metastable states, has much in common with the linear response 
formulation of transport theory within a locally harmonic approximation [^. From 
this approach one knows how the response functions have to be modified in order to 
introduce dissipation. 

2.2.1 Exploiting Green and response functions 



The time ordered average in (|20D can be identified with the two body Matsubara 
function of the one body operator F [pc| ]: 

g{q^,r-a) = -^{fF{T)F{a)),, (21) 
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On the other hand the retarded FF-response function is given by 



X 



R 



(qo, t-s)='-e{t-s) {F{t)F^s) - FHs)F{t))g, = x{qo, t - s) . (22) 



h 



(Henceforth, we shah omit the upper index "R"). It describes the response of the 
expectation value {F) to the variations of q in real time evolution, 

/CX3 
dsx{qo,t-s){qis)-qo), (23) 

The spectral representations Q{qQ,ih'r) of (p|) and xiqo^^) of (22) are obtained by 
Fourier series and Fourier transformations, respectively. As both have the same spectral 
density, one may prove [^] them to be connected by the analytic continuation 

Q{qo,ii^r) iVr ^ OJ + ie -x^iqo,^)- (24) 
The response function may be continued to the whole complex plane via H] 

P dn x"{qo,^) r T ^ n ^9^^ 
X{qo,z)= / for Imz / , (25) 



with x"iqoT^) being the imaginary (dissipative) part of x(q'O)'^)- The form ( pSD defines 
two branches. The one which is analytic in the upper half plane coincides with the 
retarded function x^iqo^z) and the one analytic in the lower half plane defines the 
advanced function x^iqo^^)- Both branches may be continued analytically into the 
other half planes. Below we will make use only of the retarded response x^iqo^^) = 
xiqoyz). On the imaginary axis {z = iw with G R) it has the following symmetry 
properties: 

ixiqo, iw))* = x{qo, (iw)*) = x{qo, -iw) = x{qo, iw) (26) 



This property, together with the relations ( ]2l| ) and (|2J) may be exploited to calculate 
the r-integrals in ([20|) as 

ln{Ug)lf^^ = P^\qrW{qo,ii^r). (27) 

r>0 

Mind that because of the reality of the collective variable one has q* = q-r- The result 
(27) may be plugged into (13) to arrive at the following form 

CP'^^(/?,go) = / V'q exp ((^^{l+kx{qo,i'^r))\qrA • (28) 

The remaining integrals hidden in Vq are of Gaussian type. As we stick to the case 
k < 0, they cause no problem as long as 

1 + kxiqo, i^r) > for r > . (29) 

As we shall see soon this leads to a condition on the temperature below which the 
PSPA breaks down, as already noticed in [12, |^]. Here, this condition only has been 
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rewritten in terms of the response functions used in the hnear response approach to 
nuclear transport (see e.g-Q). In this language the final result for C^SPA j-gg^^jg 

CPSPA(^^g^) = [](! + kx{qo,iyr))-' (30) 

r>0 

and that for the partition function of the PSPA becomes 



/ + 00 
dgo e2S'?o z{(3, qo) + kxiqo, m-))'^ ■ (31) 

-OO „^ n 



-2-nk 



2.2.2 Response functions in the independent particle model 



Before we are going to discuss further the condition ( p9| ) in the next subsection, let 
us recall how the response function looks like in the model of independent particles, 
as defined by the Hamiltonian ho of (^). It is not difficult to convince oneself of the 
following form for the dissipative part of the FF-response function 

}("{qo,uj) = -^^lF/fc(go)P riikiqo) 6{uj - eik{qo)/h) , (32) 
l,k 



where 



Fik{qo) = {liqo)\F\k{qo)) 

e/fe(go) = e/(go) - efc(5o) (33) 
nikiqo) = n{ei{qo)) - n(efc(go)) 

and n(e) being the Fermi occupation numbers 

1 + exp(p(e - ^j) 
Within this model it can easily be seen, that the xiqo: is given by 

X(?0, z) = -W \FM\' -iilf- (35) 
n'j^ eik{qo)/n-z 

Notice, please, that along the real axis the z must be chosen identical to a; + ie. For 
details about these response functions we may refer to |^. 

2.2.3 The crossover temperature 

Let us elaborate now on the convergence condition ( p9|) for the g^-integrals in 
finally to establish connection to an analogous condition which shows up when treating 



dissipative tunneling at finite temperature within the Caldeira-Leggett model [21, 22]. 
To this end the following identity is useful ||l^ 

1 + kx{qo,iVr) = , , (36) 

\Ak>A^r + Mqo)IK)^) 
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which is vahd for all r 7^ 0. The frequencies uJuiQo) appearing here are those of the 
local RPA associated to the local vibrations of the mean field around qq . They satisfy 



a secular equation |13|, which can easily be brought to the form 

l + kx{qo,z)=0 (37) 



by analytically continuing the function G{qo, i^r) to complex z by way of ( p^ ) and 
As the denominator of the ratio on the right of (|3^ ) is real and positive the condition 
( p9D can be reformulated as 

J] {ul + ujliq^)) > (38) 

V 

as already mentioned in In case that all local RPA modes are stable, and 

hence that all cjj^(go) are real, the condition is fulfilled for any temperature, viz for 
T > To = (mind (|5|)). For unstable RPA modes, on the other hand, one pair of 
corresponding frequencies w™'^*'((7o) becomes purely imaginary, in which case ( ^8|) can 
be fulfilled only above a certain minimal temperature ro(go)- The latter may vary with 
qo, but it is possible, of course, to define a minimal global temperature Tq by 

h\J-^\qo)\ 

iQ = max (39) 

27r 



such that ( [29| ) is fulfilled for all T > Tq. This temperature is identical to the so 
called "crossover temperature" (here of course for an undamped system) that shows up 
in the Caldeira-Leggett model when dealing with unstable modes of dissipative quan- 
tum systems [^]. There, the notion "crossover" indicates a transition in the nature 
of the decay of a metastable system. Above Tq the process is dominated by thermally 
activated decay ("thermal hopping") with the effects of genuine barrier penetration in 
the quantum sense to become dominant only below this Tq (called "dissipative tunnel- 
ing" for damped quantum systems). Evidently, in a typical situation, the Tq of ( p9| ) 
would correspond to that go where the top of the barrier is located. 



3 The PSPA for dissipative phenomena 

To elaborate on the connection to the treatment of dissipative tunneling within the 
Caldeira-Leggett model we need to introduce dissipation. As mentioned previously, 
the most natural way is through the response function. This can best be seen at the 
secular equation (|37|). For real z = lo = lo*, for which the response function splits into 
its real (reactive) and imaginary (dissipative) part, x{^) = x'{^) + ^x"('^)) one gets 

1 + kx'iio) = 
X (w) = . 

Whenever the function x"{^) is given by a discrete sum of (5-functions located at 
eik{qo)/fi, as shown in (^), the second equation is automatically fulfilled at the so- 
lutions L^i, of the first equation. These solutions are either real or purely imaginary 
without any sign of dissipation, reflecting the fact that the local RPA as discussed 
above corresponds to time reversible dynamics. This argument shows that irreversibil- 
ity is intimately related to the functional form of the dissipative part x"{'^) of the 
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response function. A genuinely microscopic approach would require to consider ex- 
plicitly couplings of the simple particle-hole configurations to more complicated states 



[17|. Definitely, this is beyond the scope of the functional integral method underlying 
the present model. In a more phenomenological approach one might argue to dress 
the single particle states with complex self-energies which itself may vary with tem- 
perature, for details see H] or where the inclusion of pairing is discussed. An even 
simpler way is to effectively perform the transition to a continuous spectrum, which 
directly corresponds to the procedure one employs in the Caldeira-Leggett model in 
typical solid state applications - @, [|^]. However, even for a finite nucleus such a 
transition is justified for not too small excitations. Indeed, as one knows from nuclear 
reaction theory [^] , for not too small energies resonances do overlap, implying that the 
true compound states lie dense for excitations above about 10 — 20 MeV. On the level 
of the independent particle model one simply might employ energy averages, which in 
turn are related to finite observation times of the system; for details the reader may 
be referred to (5| . 

In this paper we would not like to penetrate any further into this discussion. Rather, 
in the sequel we would like to assume the 'x"{uj) to be a continuous function of uj. In 
this case the secular equation (37) may no longer be written as in (Hol) and its solutions 



become complex quantities. To be specific, instead of (|32D we like to suggest and work 



with the following model function consisting of two Lorentzians of width r(go): 

It may be characterized as a generalization of the degenerate model often used in 
nuclear physics (see e.g. [[l5[| ) to one where the nucleonic states are spread over a certain 
region determined by the width r(go)- The strength of these intrinsic excitations is 
parameterized by the quantity -F'^(go)- A straight forward generalization could be seen 
in a summation of more than one term. In a sense the reduced form (^) corresponds to 
what has been called the "one pole approximation" (see e.g. Q). It is valid whenever 
the strength distribution is dominated by one peak, which then finally implies to have 



one prevailing collective mode. The parameters appearing in (|4l|) could be calculated 
in various ways, as indicated within the linear response approach, for instance, but 
even the Random Matrix Model (RMM) might be used (see e.g. ||5|). 

Inserting the spectral density ( |4l| ) into ( pSD the full response function can be cal- 
culated, which is needed both for the secular equation (^) as well as for the condition 



(29). The integral can be carried out with the help of the residue theorem noticing that 
the integrand has five poles altogether, situated at = z and Q, = ±£"((70) ± iV{qQ)/2, 
and closing the loop in the appropriate half plane. The final result for the retarded 
response function reads 

x{qQ,z) = F'^{qo) ° — — 2 ; (42) 

^^(^o)^ + (r(go)/2)^ - iT{qo) z- 

(The advanced response function would be obtained by changing —i into +i). For 
the condition (|29| ) one needs to know this function along the positive imaginary axis. 
There, the denominator is always positive implying that xiQa^'^w) is finite for real w. 
Furthermore, it is seen that xiQo^'^'^) still is real for continuous spectra. 
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3.1 Transport coefficients of collective motion 

We are now going to write the secular equation for collective motion in terms of trans- 
port coefficients, as it is known for the damped oscillator. This is achieved best by 
rewriting ( ^2|) in the form of the oscillator response function 

X(90, z) = / 2 . -rf I n2f \ = Xosc(go, z) . (43) 

M{qo) + iT{qo)z -il^{qo) 

The parameters introduced here correspond to the nucleonic (or "intrinsic") motion at 
any value of qo and are uniquely derived from ( ^2[ ) as follows: 

1 d^x~^{qo,z) 



M{qo) 



2 dz^ 



z=0 



£{qo) F^qo) 



(44) 



MiqM ^ x-(.o,.^0) ^ ^to); + (r(.o)/2)- 



M{qo)riqo) 



dx ^{qo,z) 



d 



z 



r(go) 



z=Q 



(46) 



These transport coefficients may be interpreted as the (local) coefficients of inertia, 
frequency and friction for the nucleonic mode. Plugging (43) into ( |37|) one obtains 



n 1^, . ^ z^+iriqo)z-nHqo)-k/Miqo) 

o = i + kx{qo,z) = 2 I -r^ \ ?^2r~\ ■ 

z^ + iT{qo)z - il^iqo) 

This equation may be fulfilled only for a vanishing numerator, which leads to the 
secular equation for the local frequencies z^(go) of collective motion, namely 

{z^)^+iT{qo)z^-w\qo)=0, (48) 

with the local collective frequency being defined as 

w^{qo) = n\qo) + k/M{qo) < n\qo) . (49) 

The last inequality is given because we are dealing with isoscalar modes where A; < 0. 
Notice that the collective frequency w may become purely imaginary, whereas the 



intrinsic one O is always real (see (45)). 

Now the frequencies 2^ (go) are no longer real quantities. A convenient form is seen 
to be: 

z^iqo) = |ro(9o)| (±Vsgnro2(go) - rf{qo) - iviqo)) , (50) 

with 

Viqo) = . (51) 

2|ro(go)| 

The dimensionless parameter r?(go) measures the degree of damping: It is smaller 
(larger) than 1 if the (local) collective motion is underdamped (overdamped). In the 
stable case tt7^(go) > the frequencies z^{qo) of (jsol) are found in the lower complex 
half plane symmetrically to the imaginary axis for r/(go) < 1 and on the negative 
imaginary axis for r/(go) > 1- In the unstable case tz7^(go) < they always lie on the 
imaginary axis, but now the frequency z~^{qo) is in the upper half plane. 

It may be worth while to briefly compare (E^ with the undamped case. It is 
easily recognized that for vanishing T the form (|47|) turns into ( ^7|) under the following 
conditions: (a) eq. ( |36| ) is evaluated at z instead of ii'r, (b) simply one (pair of) local 
collective mode(s) ro(go) is considered instead of all local RPA modes uj^{qQ), (c) the 
intrinsic frequencies e/fc(go)/^ are replaced by 0,{qo). 
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3.2 The crossover temperature for damped motion 



There is no change in the condition ( poD for convergence of the integrals in ( 28 ) . It is 
only that ( p9D takes on a different form in terms of the transport coefficients. Moreover, 
the x{QOj i^r) may be expressed by the transport coefficients by making use of (|47|). In 
this way (|2^) turns into 

i/2 + r((7o)z^r + ro^(go) >0 , (52) 



as a natural generalization of (pq). Still, for a real collective frequency this condition 
is always fulfilled. For a purely imaginary one, on the other hand, ( ^2|) can be fulfilled 
only if the Vr is larger than 



= |ro(go)| y-r]{qo) + Vv'^iQo) -sgnw'^{qo)j . (53) 
Hence, T has to be larger than the local crossover temperature 

To(go) = (x/TT^ - vi%)) ■ (54) 

Evidently, the To(go) is decreasing with growing damping strength ?/(go)- For r/(go) ^ 1 
one has ro(g'o) ^ l/2?/(9o)- The global crossover temperature, finally, has to be defined 
as 

To = maxro(go) • (55) 
For vanishing damping we recover (|39|) with w{qo) being identical to uj,y{qo). 



4 The fission rate within the PSPA 

Imagine that we are given a heavy nucleus which may decay by fission, a process 
which is to be understood as collective motion across a barrier. It is known that at 
smaller temperatures this barrier may have substructure due to shell effects. Such 
details shall be neglected here. Rather we shall assume the process to be dominated 
by just one potential minimum and one pronounced barrier. Likewise, we shall discard 
any evaporation of light particles and 7's. Moreover, the transfer of energy from the 
collective degree of freedom to the nucleonic ones will be supposed not to change much 
the latter's temperature. Under such circumstances the previously discussed path 
integral formulation may be applied, with a fixed temperature. As noted earlier, for 
the PSPA we expect great similarities to processes which are described on the basis of 
the Caldeira-Leggett Hamiltonian. 

There, the decay rate R of unstable systems at not too small temperatures is traced 
back to the imaginary part of the free energy. As can be seen in the literature, see e.g. 
1^] or [||, for T > To the following formula is in wide use 

i? = -^^Im^(/5). (56) 



It has originally been postulated in |26, ^] and, in strict sense, still lacks a general 
proof from first principles. However, it can be said that it is capable of reproducing 
correctly certain limits. For instance, one recovers correctly Kramers' high temperature 
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limit, and in the quantum limit one gets the same functional form as obtained with 
real-time path integrals Q or in a quantum transport theory Q. 

To evaluate the imaginary part of the free energy one still uses the relation J-{fi) = 
—T In Z{I3) to the partition function. For an unstable system the latter attains an 
(exponentially small) imaginary part. Following Langer |26] this may be shown by ap- 



plying the saddle point approximation and distorting the integration contour into the 
complex plane at the barrier. Expanding the logarithm to first order in the exponen- 
tially small quantity Im Z(/3)/Re Z(/3) the imaginary part of the free energy becomes: 



Plugging (^) into ( pq ) we obtain 

^"ITr^- ^^^^ 

The partition functions appearing here may be evaluated within the PSPA extending 
formula (|3^) to a dissipative system as outlined in section ^. Applying the saddle point 
approximation to the (/o-integral in (^) we obtain 

^™(/3)L = , , exp (59) 

^/-kCr{qa) 

as the contribution from the minimum and the purely imaginary expression 

ZPSPA(^)| = ' : exp(-/3.FSP^(Ag,)) CPSP^(/?,gb) (60) 

as the contribution from the barrier. Here, the stiffnesses 

d^r''''Hfi.qo)/dql=Cr{q,) (61) 

of the SPA free energy at fixed temperature appear, as it was assumed that the in- 
tegrand is dominated by the exponential and that the correction factor C^^^^(/3, go) 
varies smoothly with go- The stationary points are thus defined by this free energy 
through dJ-^^^/dqo = 0. Evaluating the intrinsic free energy in SPA from (|l9D with 
C = 1 it is easy to convince oneself that the extremal points fulfill the relation 

qa/b = k{F),^^^ . (62) 

The derivatives of the eigenvalues €i{qo) with respect to go needed here may be obtained 
from time-independent perturbation theory. In (|6^ ) the indices a and b stand for the 
minimum and the maximum (or barrier) of .F^^^(go), respectively. The relations ( |6^ ) 
are nothing else but the self consistency condition (|l^) applied to the two stationary 
points of the system. Whereas (|5^ ) was obtained through the common Gaussian in- 
tegrals of the saddle point approximation, for (l60| ) the integration contour had to be 
deformed such that it runs parallel to the positive imaginary axis. This is the reason 



for the additional factor 2 in the denominator of (|60|) , see |27|, ^] . 

The generalization of the PSPA correction factor (^) to damped quantum systems 
may be performed by replacing the response function (|35D of the independent particle 
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model by its continuous version ([4 2D. Furthermore, we may make use of the transport 
coefficients introduced in section |3.1| . In this way one gets: 

,2 



*) = n -;+^'^°'''- + . (63) 



As mentioned earUer, in comparison to (|36| ) there is only one (pair of) mode(s). The 
relation between the local nucleonic frequency ^{qo) and the local collective frequency 
^(go) is given by (^9|) . It is worth stressing that the infinite product (|63|) is convergent. 
To guarantee this important feature, it suffices to have the same coefficients for local 



inertia and damping in the numerator and the denominator [29|. 

Plugging (|4|), (^9|) and (|6^) into (|58| ) we obtain the following expression for the 
PSPA decay rate of the system under consideration: 



Like in the sequel we have partly used indices "b" instead of an argument qb to keep 
our notation short. The two first factors, which in a sense represent dynamics, have 
come in through the crossover temperature ro(go) discussed in section p.2| , mind ( |5^ ) in 
particular. Like in the Caldeira-Leggett model the decay rate factorizes into a classical 
part R^^^ and a quantum correction factor f^^^^ 



R ^ — -^class ^ /qm ) (65) 



with 



and 



jPSPA _ C^^^^{P,qh 
respectively. 

Let us discuss first the factor -R^iass^ which survives the classical limit, as no h is 
involved. Evidently, it contains the common Arrhenius factor 

exp (-/3Sb) = exp (-/3(^SPA(^^ _ _^SPA(^^ ^^^^^ ^ ^gg^ 

defined here by the difference £'b of the free energy between barrier and potential 
minimum. The influence of damping is given by the correction factor found first by 
Kramers l3^, namely 



fK = ^Jl + vt-Vb . (69) 

It decreases monotonically with increasing % and for % ^ 1 behaves like 1/2%. The 
remaining factor can be made to become proportional to the attempt frequency Wa at 
the minimum by writing 



27r V \CT{qb)\ 27r V Mfe V C^'^'Kqa) V 
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Here, use has been made of the relation between the frequency and inertia of the local 
mode and the associated stiffness: 



C'^°'\qo)=M{q^)w\qo) (71) 
Putting all factors together the classical rate may be written as 

It contains Kramers' original form 



RK = ^e-P''- (VlT^-^.) (73) 
as the first factor. In addition there is a another correction factor, 



.PSPA^ K l \C^°\%)\ / CAqg) 

^^1-^ V V C™"(g„) V \C^{q^)\ ^ ^ 

not present in the derivations based on the Caldeira-Leggett model. The reasons for 
that are obvious. First of all, in this model the inertia in the collective mode is simply 
put equal to a constant which renders the first factor on the right of (^4|) equal to 
unity. Second, the dynamical stiffness (|7l| ) is forced to be identical to the one of the 
phenomenologically introduced collective potential, d'^V{q)/dq^ in [||, ^, ^. This is 
achieved by working with a Hamiltonian in which from the beginning the collective 
part is renormalized by the term x(0)(7^/2 which in the linear response approach is 
induced by the static influence of the coupling. In this way only the dynamical part 
of the induced force appears, which in the end may lead to Ohmic friction. In our 
approach, where all transport properties of the collective dynamics are generated from 
the two body interaction, such manipulations are not meaningful. In certain limits it is 
possible, however, to simplify the f^-^J^ of (74). For slow collective motion, sometimes 



referred to as the zero frequency limit, the local stiffness ( |7l|) of collective motion may 
be shown to be represented by that of the free energy (^Tj) (see e.g.||5|) such that one 
simply has 



fciZ^ ~ y ^ • (75) 

For a derivation of this factor based on a generalized version of Kramers' equation and 



picture of the decay we like to refer to |31] 



Next we turn to the quantum corrections to the classical rate (pq), which in this 
approach is given by the ratio of the PSPA corrections C^^^'^(/3, go) of (|63|) evaluated 
at the barrier and the minimum. With the help of the local nucleonic and collective 
frequencies r^(go) and voi^qo) it writes 

j^PSPA _ TT + Tfe Vr + ^1 Y\ + ^r- + 

The nice feature about this structure is that it converges for all conceivable values 
of the transport coefficients as long as T > Tq. The reason simply is that it is the 
ratio of two convergent products of type (|63|). As the alert reader may guess a simple 
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generalization of the quantum correction factor of the Caldeira-Leggett model (see e.g. 
^]) to coordinate-dependent coefficients like 

.LC , TT + T^r + , . 

(see e.g. [^2|) may (for Ohmic damping where T does not fall off for large frequencies) 
lead to problems of convergence. Indeed, for Ta 7^ the infinite product either 
converges to zero or diverges depending which one of the F's is larger p^|. In [H] the 



form (77) has been derived on the basis of a quantum transport equation. In |33] 



this factor has been evaluated for microscopically calculated transport coefficients. To 



circumvent the convergence problem in (77) the individual F's had been replaced by 
the arithmetic mean value 2? = + r^. 

The local frequency of the collective motion w is real at the minimum and purely 
imaginary at the barrier, whereas the frequency of the nucleonic motion is real 
everywhere. The denominator of the first term in ( [7^ ) vanishes as T approaches Tq, 
corresponding to the definition of the crossover temperature in section ^ (see (^2|)), 
but all other factors are strictly positive. For this reason, at Tq the quantum correction 
factor /™^'^ diverges to plus infinity, a feature well known from the Caldeira-Leggett 
model p^ . 

In the limit of very high temperatures T 3> /iJl, where h£t represents the typical 
nucleonic excitation, /™^^ strictly converges to unity. For nuclear fission collective 



motion is expected to be slow in the sense [31| of having 



41MeV 

1 MeV « /itu < no -— ^ 6MeV (78) 

One may expect quantum effects in collective motion to disappear already for T ^ hvj. 
Indeed, the quantum correction factor will be close to unity already for T/Kl = 0(1). 
At least this can be shown for the two factors of ([7^). Divide all numerators and 
denominators by v"^ and neglect {vo/vr)'^ <C 1 in the resulting denominators. This is 
justified simply because the condition T/Ml ~ 0{\) implies O/f,. ~ l/(27rr) (mind 
(H)). Remain the terms which involve F/f^- They can be neglected if T ^ l/27r • ^F 
or T/hw ^ ^/tt. Microscopic computations of the transport coefficients show this 



condition to be fulfilled, although r] itself increases with T; see figure 3 of fSlI] or figure 
5.2.10 of llEtl. 



5 Conclusion 

We have been able to demonstrate how the PSPA can be extended to treat the decay of 
damped meta-stable systems. In this first step a simple schematic two body interaction 
has been taken and the nucleonic excitations have been assumed to be concentrated 
in one Lorentzian peak around a certain mean value. Generalizations to more general 
systems should not cause too many problems. So far we concentrated on the quantum 
corrections to thermal hopping which take place above the critical temperature Tq 
[p2|, |6|. At this temperature the common semi-classical treatment of functional integrals 
breaks down, simply because for unstable modes the Gaussian integrals diverge for 
smaller temperatures. So far this latter feature also limited the applications of the 
PSPA to bound systems |]l^, 12, 13|. There is hope that this deficiency can be overcome 
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in very much the same way as it was possible for dissipative tunneUng 34]. Work 
in this direction is under way ||35[| . 

There are several advantages of the method presented here, both over the usual 
approach to dissipative tunneling within the Cladeira-Leggett model [Q, as well as 
with respect to the Locally Harmonic Apprtoximation (LHA) |^] to quantum transport. 
Different to the Caldeira-Leggett model, all transport properties derive from the two 
body interaction of the many body system. No phenomenological assumptions have 
to be made for any transport coefficient. The effects of the two body interaction are 
treated on a fully self-consistent level, largely because the collective variables can be 
introduced globally by way of the Hubbard-Stratonovich transformation. For the LHA, 
on the other hand, and on a quantum level this is possible only locally ||5[. This method, 
however, is more flexible with respect to the thermal properties. There, one needs not 
rely on the concept of a fixed temperature, an assumption which is questionable for 
isolated systems. 
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